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Abstract 

Graphene is a promising candidate for optoelectronic applications such as photodetectors, 
terahertz imagers, and plasmonic devices. The origin of photoresponse in graphene 
junctions has been studied extensively and is attributed to either thermoelectric or 
photovoltaic effects. In addition, hot carrier transport and carrier multiplication are 
thought to play an important role. Here we report the intrinsic photoresponse in biased 
but otherwise homogeneous graphene. In this classic photoconductivity experiment, the 
thermoelectric effects are insignificant. Instead, the photovoltaic and a photo-induced 
bolometric effect dominate the photoresponse due to hot photocarrier generation and 
subsequent lattice heating through electron-phonon cooling channels respectively. The 
measured photocurrent displays polarity reversal as it alternates between these two 
mechanisms in a backgate voltage sweep. Our analysis jdelds elevated electron and 
phonon temperatures, with the former an order higher than the latter, confirming that hot 
electrons drive the photovoltaic response of homogeneous graphene near the Dirac point. 



Optoelectronic properties of graphene have attracted substantial interest due to 
graphene's high carrier mobility, zero bandgap, and electron-hole symmetry/ Graphene 
can absorb light and turn it into a photocurrent over a wide range of the electromagnetic 
spectrum, from the UV to the visible and infrared regimes. A photocurrent response at 
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graphene-metal contacts was reported by various groups including ours, " and in these 
early reports the origin of the photocurrent generation was attributed to the classical 
photovoltaic effect. The thermoelectric Seebeck effect, induced through photo-excited 
carriers, was tirst identitied in single layer - bilayer graphene junctions experiments as 
the origin of photocurrent.^ The technologically important photocurrent of graphene p-n 
junctions has been either described as photovoltaic ' or thermoelectric " in nature. Both 
mechanisms may be enhanced by hot carriers^ ^"^^ that persist after photoexcitation in 
graphene due to an elecfron-acoustic phonon decay bottleneck.^ Elecfron-elecfron 
scattering on the other hand dominates the photocarrier energy relaxation, where the 
energetic photoexcited carriers multiply while decaying toward an elevated electronic 
temperature Te.^ '^ The decoupling of electron and lattice temperatures in graphene may 
also explain how thermoelectric confributions can be compatible with the fast 
photoresponse on the order of lOOGHz seen in modulated photocurrent 
measurements. ' Due to the identical polarity of photovoltaic and thermoelectric 
currents in metal-graphene or graphene p-n junctions, it has been quite a challenge to 
identify the dominant photocurrent mechanism. However, the two effects show opposite 
polarity in the unipolar junction regimes p^p and n"n, which allowed the identitication of 
the thermoelectric effect in the photocurrent of graphene on atomically flat boron-nitride 
dielectrics."^ Very recently, a bolometric photoresponse was observed in gapped bilayer 
graphene.'^ 

In this work, we investigate the origins of the photocurrent response in graphene 
by measuring the photoconductivity of the homogeneous graphene channel in a graphene 
field-effect transistor (FET) fabricated on Si/Si02. In this simplest and yet unstudied 
configuration, the photocurrent polarities due to photovoltaic and thermoelectric effects 
are opposite, which allows us to directly determine the dominant photocurrent generation 
mechanism. In addition, the thermoelecfric effect is an order of magnitude smaller than 
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the photovoltaic effect observed experimentally. Our results clearly reveal the importance 
of the photovoltaic effect particularly near the Dirac point, while a photo-induced 
bolometric effect dominates in the n-type and p-type doping regimes. 

The graphene FETs in this work are biased at one contact with a moderate drain 
voltage on the order of = -IV . The doping level in the graphene is controlled 

electrostatically via a global silicon backgate. While a chopped and focused laser beam at 
wavelength 690nm is scanned over the sample, we measure the AC photocurrent 
amplitude and phase, which are obtained by a lock-in amplifier referenced to the 
chopping fi-equency (Fig. la). Photocurrent images (Fig. Ic) show a uniform photocurrent 
within the interior of the graphene channel, exhibiting a photoresponse on the order of 
10 A/ W . Contact effects, which have been discussed extensively for the zero bias case 
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(Fig. lb), " are limited to the immediate vicinity of the leads. Under applied bias, a small 
photocurrent (~10% of the observed value at graphene) is also observed when the laser 
beam is positioned a few microns away fi-om the graphene, over the bare Si/Si02 
substrate. The photocurrents from graphene and substrate exhibit opposite signs at 
Fq = 5V , as can be seen in the phase image. We show in the Suppl. Info, that the 

substrate-related photocurrent is due to a 2mV photovoltage generated at the Si/SiOa 
interface that is sensed by the graphene field-effect transistor. This "photo field-effect" 
switches sign when the Fermi level in the graphene is swept past the charge neutrality 
point. We can easily correct for its effect at the graphene position. 

When photons are incident directly on the graphene channel, a true photocurrent 
from graphene ensues. This photocurrent switches sign twice when the Fermi level in the 
device is swept fi-om p-type over to intrinsic and then to n-t5^e (Fig. 2a). Figures 2b and 
2c show the photocurrent in the center of the graphene as a fimction of gate voltage. 
Between gate voltages of Vg=OV and Vg=2.3V, which is close to the Dirac point, the 
photocurrent is positive and points in the direction of the DC current (from source to 
drain), while outside of this range, the photocurrent is negative. The largest photocurrent 
amplitude is reached at high electrostatic doping (positive or negative). Figure 3 shows 
the full gate and drain- voltage dependence of the AC photocurrent near the center of the 
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device. The six-fold pattern observed here is reminiscent of the situation observed in a 
dual-gated graphene p-n junction without source-drain bias, where a six-fold pattern as a 
function of the two independent gate voltages has been attributed to the thermoelectric 
Seebeck-effect.^'^*^ In graphene under bias, the six fold pattern is composed of a two-fold 
sign change with gate voltage, and a single sign change with drain voltage. The well 
known tilt in the current-voltage characteristic (Fig. 3a) due to channel doping by the 
drain voltage, is also apparent in amplitude and phase of the AC photocurrent (Fig. 3b, c). 
The drain voltage dependence of the photocurrent is linear up to voltages on the order of 
a Volt, above which saturation sets in, especially at high doping (Fig. 3d and Suppl. 
Info.). 

Several processes of photocurrent generation may play a role in biased graphene: 
(1) The thermoelectric Seebeck effect (TE): The laser spot produces a temperature 

gradient in the device, which together with a doping asymmetry generates a 

thermoelectric current. A typical experimental setup involves a p-n junction to produce 

Sin 

the Seebeck effect. ' While unbiased homogeneous graphene is not expected to produce 
a thermoelectric current, the presence of an applied drain voltage renders the 
electrochemical potential spatially variable, and a TE effect should ensue even in 
homogeneous graphene; (2) The photovoltaic effect (PV): Photo-excited electrons and 

holes are accelerated in opposite directions by the electric field. They produce a 
photocurrent either by reaching the contacts still hot, or by establishing a local 
photovoltage within the focal area that drives the photocurrent through the rest of the 

Q 

device; (3) The bolometric effect (BOL): The incident electromagnetic radiation raises 
the local temperature of the graphene, which alters the resistance of the device, producing 
a change in DC current under bias. Fig. 4 (a-c) summarizes the relative magnitude and 
signs of these effects in non-biased and biased graphene photodetectors. In the following 
we will discuss these effects in greater detail. 

It is instructive to first review the basic TE effect in a p-n junction at zero 
bias.^'^'^*^ Figure 4a illustrates the typical experimental situation, where the laser spot 
induces a temperature gradient and results in a net thermoelectric voltage across the 
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graphene due to the Seebeck effect. The generated photocurrent is proportional 
to (/Sj - /Sj ) dTj dx , where is the Seebeck coefficient for the two sides of the junction. 

In the zero bias case, the sign of the thermoelectric photocurrent is the same as the 
photovoltaic effect, except for some unipolar junction regimes. This feature is most 
pronounced in very clean samples, leading to the six-fold sign change recently observed 
in graphene on boron nitride substrates.'*' When the graphene doping is homogeneous, the 
photocurrent will be zero. However, the uniform doping can be rendered asymmetric 
under an applied drain bias, such that the effective doping along the graphene channel 
changes gradually. This doping asjmimetry and the associated Seebeck effect are more 
prominent at low doping. Typical energy band profiles are illustrated in Figs. 4b,c 
together with the sign of the various photocurrents relevant in our experiments. Near the 
charge neutrality point, the TE registers a photocurrent opposite in sign to both the 
experimentally measured photocurrent and the expected photovoltaic current. Since the 
Seebeck coefficient is related to the electrical conductivity cr through the Mott formula, 

7T^ k^T 1 d iJ ^ 

S = , application of the thermoelectric theory yields us an estimate of a TE 

3e <j ds 

photocurrent of approximately 4nA near the Dirac point (see Suppl. Info.), opposite and 
an order smaller than that observed in experiment. Away from FG=fDirac, the TE effect is 
expected to be even less important because the electrochemical potential asymmetry is 
reduced with increasing Fermi level. 



Next we consider the bolometric response of the graphene photodetector, which 
can only be observed in biased devices, and was therefore absent in previous short-circuit 
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measurements on graphene contacts or junctions. " The bolometric effect can be 
determined by a measurement of the temperature dependence of the transport current. 
Figure 5a shows the change in transport current M compared to its room temperature 
value when cooling down in the cryostat. From this we can extract the gate-voltage 
dependent bolometric coefficient y5(FQ) = A/(FQ)/Ar(FQ) (Fig. 5b, red curve), which 

describes the sensitivity of the transport current to changes in temperature around room 
temperature.'^"^' A similar measurement around our lowest temperature of ro=205K 
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yields the other (blue) curve for the temperature coefficient in Fig. 5b. The negative 
values of ^{Vq) in single-layer graphene are due to electron-acoustic phonon scattering 

and electron-remote phonon scattering of surface polar phonons in the SiOa underlayer, 
which both are enhanced at elevated temperatures.'^ For comparison, Fig. 5c shows the 
experimental photocurrents at room temperature and J'o=205K. The similar shapes of the 
curves in Figs. 5b and 5c are striking. In particular, both the bolometric coefficient and 
the photocurrent are largest (and negative) at high doping and the electron-hole 
asjmimetry of the temperature coefficient is reproduced in the photocurrent measurement. 
While the bolometric effect can explain these general features of the gate-voltage 
dependence of the photocurrent, it caimot explain the positive photocurrent in the 
direction of the DC current under low doping. 

The photovoltaic effect (PV) is the only known mechanism that is consistent with 
the positive photoresponse near the Dirac point seen experimentally. Under light 
excitation, the hot carriers with separated electron and hole chemical potentials are 
formed on a timescale of the order of 100 fs, followed by a slow pico-second carrier 

recombination and cooling.''"'^' At steady state, the photo-induced carrier density n^^j^ 

depends on its non-equilibrium carrier temperature and the chemical potential and has 
to be determined by imposing charge conservation i.e. nl=nl. In addition, has its 
maximum at the focal spot and is proportional to QL/k^ where Q is the laser power and 
is the electronic thermal conductivity, related to the electronic conductivity through 
the Wiedemann-Franz relation. Making use of the above facts allows us to compute the 
gate dependent photo-induced carrier densities «*/h('^G) (Suppl. Info). The photovoltaic 

current can then be estimated fi-om qn /j.^ , where ^ is the electric field and 

// « 2700cm^/Vs our device's carrier mobility. In general, the non-equilibrium 

decrease with increasing equilibrium doping, consistent with the fact that electron- 
electron scattering which serves as efficient means towards carrier equilibration increases 
with doping. Hence, the modeled photovoltaic current as shovm in Fig. 5c decreases 
gradually with increasing bias away fi-om the charge neutrality point. Near the Dirac 
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point, the elevated electron temperatures are T^=T^-Tf^ x8K and 12K at initial 
temperatures = 300K and 200K respectively as shown in Fig. 5d, which also decrease 
with increasing equilibrium doping. 

Next, we subtract the estimated photovoltaic current component from our 
experimentally measured photocurrent, to obtain the purely bolometric contribution 1^^^ . 

The elevated lattice temperature =^ph~^o then be determined from the ratio 

I^Qi^lfi, where is the experimentally determined bolometric coefficient (Fig. 5b) 

discussed earlier. As shown in Fig. 5d, is of the order of a Kelvin and is weakly 

dependent on doping. Both and are enhanced at lower temperatures due to less 

efficient cooling. Although electron-electron interactions result in an energy equilibration 
of the electronic system, they do not lead to a net energy loss. The dominant energy loss 
pathways are due to phonons. Intrinsic energy loss channels due to acoustic/optical 
phonons are known to be inefficient energy loss pathways and have only a power density 
2; ^lO'-lO" W/m' for 7^-7^=1 OK (Suppl. Info.).^^'^^ With a typical graphene- 
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substrate thermal resistance of r^^ «10" Km /W, ' our measured lattice temperature 
suggests a lattice heating power density of the order of = ^L/r^^ « 10^ W/m^ or 

somewhat lower, because of the uncertainty in the temperature at the Si/SiOa interface 
after photoexcitation. The inferred larger power density in experiments may indicate the 
presence of more efficient energy loss pathways in graphene devices. Electron-remote 
phonon scattering of surface polar phonons in the Si02 underlayer, ' " and extrinsically 
enhanced phonon-scattering processes mediated by ionized impurity and random strain 
induced gauge fields^^ may provide these additional decay pathways. More systematic 
studies in this direction are needed to clearly identify the energy dissipation mechanisms 
in graphene. 

In conclusion, our measurement and theory of the intrinsic photoconductivity of 
biased graphene show that both bolometric and photovoltaic effects are relevant in its 
photoresponse. The polarity and magnitude of the photocurrent can be modulated by 



7 



electrostatic doping, which upon examination allows us to probe the non-equilibrium hot 
carrier and phonon characteristics. Our studies therefore open up the possibility of 
engineering the hot carrier photoresponse, which plays an essential role in applications 
such as bolometers, calorimeters and photodetectors. 

Methods 

Graphene field-effect transistors are fabricated by mechanical exfoliation from 
graphite. E-beam lithography and evaporation of 15nm/70nm of Ti/Au are used for the 
contacts. A second e-beam lithography step and plasma etching defines rectangular 
graphene sheets with length and width of 6|jm and Ifjm. Single-layer graphene is 
identified by its optical contrast on the 90nm Si02 dielectric. We use a highly doped (n- 
t5^e) silicon substrate to avoid an excessive photo field-effect, which dominates 
measurements on samples fabricated on lightly doped silicon. The as-prepared samples 
behave p-type. To render them intrinsic, enhance mobility, and reduce hysteresis, the 
samples are evacuated for a few days prior to the measurements and subjected to a 
combination of current annealing and laser scanning with a green laser for an hour. After 
this treatment, the Dirac point is near =0V . AC photoconductivity measurements are 

performed with a Ti/Sapphire infrared laser at around X = 690nm . The laser spot 
diameter on the sample is about d « 700nm . Incident laser power values are given in the 
Figure captions. Considering the elecfric-field enhancement due to the gate stack, we 
estimate that 2.5% of the incident light is absorbed in the graphene. The photocurrent 
amplitudes reported in this paper are peak-to-peak values throughout. We confirm that 
the chopping frequency / = 1 . IKHz is not too large by comparing our AC photocurrent 
results with equivalent DC measurements, which have much worse signal/noise, but 
agree with the photocurrent phase reported in this paper. The photocurrent measurements 
are performed in an optical dewar with vacuum on the order of 10"^- 10"^ Torr, eliminating 
the effects of photodesorption of oxygen, and allowing us to measure the temperature- 
dependence. 
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Figure 1 




Figure 1: AC photoconductivity of graphene with and without bias, (a) Schematic of 
the photoconductivity measurement setup with /I = 690nm laser wavelength, 
P = 220|a.W laser power, d » VOOnm spot diameter, and / = 1 . IKHz chopping 

frequency, (b) AC photocurrent at Vj^=OV drain bias (short-circuit photocurrent) and 

FQ=5Vgate bias (n-type regime). The laser-scanning image (left) shows source and drain 

electrodes. The outline of the graphene sheet is indicated. Middle: Amplitude (R) of the 
photocurrent. Right: Phase ((p) of the photocurrent. A localized photocurrent is generated 
close to the contacts. The phase at source and drain contacts are (p « 180° (corresponding 
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to negative /ph) and » 0° (positive /ph) respectively, (c) AC photocurrent under a drain 
bias of Fd=-1V and gate voltage Fq=5V. Laser scanning, amplitude, and phase images 

are shown. The entire graphene channel is producing a photocurrent under applied bias, 
and the photocurrent is uniform in both amplitude (/?«46nA) and phase (^3«180°) 
throughout the channel except for the contact regions. The small photocurrent generated 
next to the graphene channel is in response to a photovoltage at the Si/Si02 interface that 
couples capacitively to the graphene field-effect transistor. 
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Figure 2: AC Photocurrent as a function of gate voltage, (a) Spatial images of the 
photocurrent amplitude and phase as a function of gate voltage. (Drain voltage 
= -IV ; laser power P = 370fa,W). The photocurrent in the graphene channel switches 

sign twice from negative {(p»\SO°) to positive {g)»0°) and back to negative 
(^wl80°). (b) Photocurrent amplitude (red) and phase (blue) in the center of the 
graphene channel as a function of gate voltage. The sign of the photocurrent is indicated. 
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For comparison, the DC current has positive sign for = -IV . (c) Photocurrent (red) 

corrected for the photo field-effect (Suppl. Info.) and source-drain current (black) as a 
function of gate voltage. 
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Figure 3 




Gate (V) Drain Voltage (V) 



Figure 3: Gate and drain voltage dependence of the photocurrent. (a) DC current, (b) 
photocurrent amplitude, and (c) photocurrent phase, as a function of gate and drain 
vohages. SoUd Unes indicate the modeled sign change, (d) Gate voltage characteristics of 
the photocurrent for different drain vohages, corrected for the photo field-effect. Solid 
Unes are Lorentzian fits to the data. A 16nA photocurrent near the Dirac point at Fd=OV 
is due to residual doping, (e) Peak bias-induced photocurrents at high and low doping as a 
function of drain voltage. (Solid lines are a guide to the eye.) The photocurrent at high 
doping is due to the Bolometric effect, while the photocurrent at low doping is due to the 
Photovoltaic effect. 
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Figure 4 




Figure 4: Thermoelectric, Bolometric, and Photovoltaic components of the graphene 
photocurrent, (a) Energy band diagram of a graphene p-n junction used in typical 
photocurrent experiments, with laser incident at the junction interface. The DC current 
and thus the BOL current are zero. TE and PV components point in the same direction, 
(b, c) Energy band diagrams for our experiments, where a uniform graphene channel is 
subjected to a finite source-drain bias Fd, under small and large electrostatic doping 
conditions respectively. The directions of the photocurrent components are indicated. At 
low electrostatic doping the PV effect dominates, while at high electrostatic doping the 
BOL effect overshadows both PV and TE effects. 
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Figure 5: Magnitude of BOL and PV effects and determination of the hot electron 
and lattice temperatures, (a) Temperature-dependent change in current from its room 
temperature (7'ti=299K) value, as a flinction of gate voltage. Inset: Gate-voltage 
characteristic as a function of temperature when cooling down from room temperature, 
(b) Bolometric coefficient /3 {V^^^ = ISJ {Vq) j l^T {y^) of the transport current as a 

function of gate vohage at two different ambient temperatures To as indicated, (c) 
Experimental photocurrent (corrected for the photo field-effect) and modeled 
photovoltage component of the photocurrent. The PV component is important only near 
charge neutrality where it determines the sign of the overall photocurrent. Away from the 
Dirac point, the bolometric effect (b) dominates, (d) Electron {T^ and lattice (7l) 
temperatures extracted Irom the data as detailed in the text. 
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I. EXTRACTING AND CORRECTING FOR THE PHOTO FIELD-EFFECT 



We are interested in the photocurrent that is generated by photons absorbed in the 
active channel of the graphene photodetector. These photons produce electron-hole pairs 
in the graphene, which rapidly decay into a cloud of hot electrons and holes, leading to 
photocurrents due to the photovoltaic, thermoelectric, and bolometric effects. In addition, 
there exists a photocurrent contribution that is extrinsic to the graphene photodetector, 
and which we would like to correct for. This contribution is due to light absorbed in the 
Silicon substrate close to the Si/Si02 interface, producing a photovoltage at the interface, 
which is picked up by the gate-sensitive graphene field-effect transistor as a change in 
source-drain current. It should be possible to avoid this "photo field-effect", by using 
metallic gates, but as we show below, it is also easy to correct for the effect because the 
intrinsic and extrinsic photocurrent contributions can be spatially decomposed. 

Due to a workfunction mismatch between Silicon and Silica, the conduction and valence 
bands in Silicon bend at the interface. For n-type doping of the Silicon substrate as in our 
case, the bands in Silicon bend upward, which leads to a triangular potential well for holes 
at the interface [T]. Photo-generated holes diffuse toward the interface, while electrons are 
repelled from the interface. This leads to an additional positive voltage on the interface, 
which acts just like an applied positive gate voltage would in the graphene field-effect 
transistor, altering the source-drain current. Since the transconductance of a graphene 
field-effect transistor switches sign at the Dirac point, the photo field-effect also switches 
sign at the Dirac point {Vcnp ~ IV in Fig. Sla). This is in contrast to the intrinsic 
photocurrent, which switches sign twice, as discussed in the main text. 

The magnitude and spatial extend of the photo field-effect depends on the substrate 
chemical doping. For intrinsic or lightly doped silicon, the carrier lifetime is long, and the 
magnitude and spatial extend can be large (centimeters). For heavily doped Silicon, as in 
our case, the lifetime is shorter, but we still measure a photo field-effect, as can be seen 
from Fig. Sla, where the photocurrent is plotted as a function of gate voltage and position 
perpendicular to the graphene channel. The intrinsic photocurrent components decay 
rapidly once the laser spot moves away from the graphene, but the photo field-effect remains 
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up to a distance of several microns. This behavior allows us to estimate the magnitude of 
the photo field-effect at the position of the graphene by considering the photocurrent that 
is generated away from the graphene and fitting it spatially to Lorentzians as exemplified 
in Fig. Sib. Figure Sle shows the values of the extracted photo field-effect at the center of 
the graphene as a function of gate voltage. As expected, the curve is proportional to the 
transconductance Qm extracted from the I — Vq characteristic. The proportionality factor is 
2nA///S at a laser power of 370 //W. This means that a photovoltage of 2mV is generated 
at the Si/Si02 interface. We can now subtract the photo field-effect component from the 
total photocurrent and obtain the intrinsic photocurrent in Figs. Sic and Slf. This latter 
result is used as the basis for our model on the photovoltaic and bolometric components of 
the intrinsic photocurrent. 

II. SPATIAL DISTRIBUTION OF THE AC PHOTOCURRENT AND PHO- 
TOCURRENT SATURATION AT HIGH BIAS 

The spatial distribution of the photocurrent in biased graphene along the channel direc- 
tion is shown in Fig. S2 as a function of gate voltage for different drain voltages. At zero 
drain voltage, the well-known contact effect is present, where regions close to the metallic 
leads become photoactive because of band-bending there. Both the photovoltaic effect and 
the Seebeck effect hkely play a role in this regime. The contact effect is strongest with the 
graphene channel electrostatic doping opposite to the metal-induced doping of the graphene 
beneath the leads, which produces two back-to-back p-n junctions. In our case the metal 
dopes the graphene n-type and p-n junctions exist for negative gate voltages. These junc- 
tions move further into the channel for gate voltages that approach the fiat-band voltage at 
Vg—2Y. At more positive gate voltages, no p-n junctions exist, and the photocurrent from 
the contact regions is smaller and is generated right at the contacts. 

Once a drain bias in excess of about Vd^O.SV is applied, the bias-induced photocurrent, 
which is the topic of this paper, dominates. The high spatial uniformity of this photocurrent 
is apparent at Vc)=-lV, where the middle 4/im of the 6fim long graphene shows essentially 
the same photocurrent and gate- voltage dependence. Contact effects are limited to a l^m 
area next to the metal leads. There is a slight tilt in the gate-voltage characteristic due to 
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drain-voltage induced doping of the channel interior, which affects the right (drain) side of 
the device more than on the left, and which shifts the photocurrent pattern down by IV 
at the drain and half of that (0.5V) in the center of the device. This tilt becomes stronger 
at Vd=-2V and -3V as expected. In fact, one can use these photocurrent measurements to 
determine the Dirac point inside the biased graphene channel as a function of x-position. 

The saturating behavior of the bolometric component of the photocurrent is already 
becoming apparent below Vd^-IV (see main text Fig. 3e). The color-scale bars in Fig. S2 
show that at higher drain voltages both the BOL and PV components indeed saturate. Once 
the electron temperature is elevated due to the bias, additional photogenerated carriers will 
not be able to increase the electron temperature as much as before, because the photocarrier 
lifetime will be reduced if the electron distribution is already hot. The high bias thus limits 
both BOL and PV components of the AC photocurrent. 



III. DEVICE MODELING 

We consider back-gated (Vg) graphene devices, where the left contact is grounded i.e. 
Vl = and Vr allowed to vary. Our model considers the operating regime where the bias 
current Idc induced by Vr is still in the linear regime. The electrochemical potential jj, in 
the graphene channel {—^<x<^) is simply, 

l^i^) = - ^ (1) 

The electrical potential energy $(x) (or Dirac point energy) is given by, 

$(x) = ^^a; + ^^+/.(x) 

Pl/r = -sign(VG - Vl/r) X hvf^lnCB \Vg - Vl/r\ (2) 

To keep the analytics tractable, we fit the electrical conductivity phenomenologically for 
electron-hole puddles, 

(7(e) = ^v/?+A^ (3) 

where e is defined to he e—ji — ^. amin is the minimum conductivity and A represents the 
neutrality region energy width. Both can be simply extracted from the experiments, through 
(a) = J a{e)dx. In our experiments, device physical dimensions are W x L — 1 x 6/im 
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and tox = 90 nm. The experimentally measured graphene electrical conductivity is fitted to 
Eq.jsj with best fit values of (Jmin = 2.3 x 10^^ S and A = 75meV. In our experiment, the 
extracted effective mobility around the neutrality point is /z = 0.27m^/Vs. 



IV. THERMOELECTRIC CURRENT MODELING 



The Seebeck coefficient is computed using the Mott formula[2], 

TT^klT 1 da vr^HT 2e^ 



' 3e ade 3e + ^ ' 



The second equality makes use of Eq.|3} Hence, Sg for each location in x can be computed. 
The photocurrent density (Am~^) generated by the thermoelectric effect can be computed 
through, 

L 

Jte = -^ ,Sg{x)—dx (5) 

As mentioned in the main manuscript, the uniform channel doping can be rendered 
asymmetric under an applied drain bias, such that the effective doping along graphene 
changes gradually across the two contacts. This spatial variation in doping is described by 
— (see Eq. [5]), from which the resulting Seebeck coefficient can be computed from 
Eq.i 

Consider photo-excitation in the middle of the graphene channel. A simplified model for 
the hot electron/hole temperature profiles due to photo-excitation suffice [3j: 

re/h{x) = ^k{x) + To (6) 

where To is the ambient temperature, Q is the absorbed laser power, L the device length and 
kq is the electronic thermal conductivity, where kq and a are related through the Wiedemann- 
Franz relation. A(x) is a triangular function, with maximum at the middle of the channel i.e. 
X = and zero at x = ^\L. As discussed in Sec. |v| Te/h{x) has a maximum temperature of 
8 K in the middle of the channel. The thermoelectric current calculated from Eq. |5] yields 
Ite ^ 4nA at Vr = 1 V and when graphene channel is biased near charge neutrality. This 
thermoelectric effect is an order smaller than the corresponding photocurrent observed in 
experiment and also has an opposite sign. 
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V. PHOTOVOLTAIC CURRENT MODELING 



As argued in the main manuscript, the observed photocurrent of Ipy ~ 40 nA in graphene 
when biased near the charge neutrahty point (i.e. Vq = 0) is due to a photovohaic contri- 
bution. The photovoltaic current can be modeled by, 

Jpv = cT*C =^Wr-Pl- eVji) (7) 

where a* is the photoexcited conductivity. With an applied drain bias of Vr = —IV and 
source Vl = V, the calculated channel electric field (using Eq. [2]) when the device is biased 
near the charge neutrality point is ^ = 1.53 x lO^V/m. This yields us a* = 2.6 x 10^^ S. 
Since a* can be expressed as a* ~ qn*fi*, where n* is the photo-induced carrier density and 
H* the effective mobility of these excited carriers, where we assumed fi* ^ fi = 0.27m^/Vs, 
where fi is inferred from experiments. We obtain n* = 6 x 10^^ m~^ at = 0. 



The photo-induced electron and hole densities at the laser spot are estimated to be 



n' = n 



*^ ~ QfTi* /2, where Qf = WL/aspot ~ 16 is a geometrical scaling factor with agpot 
being the focal area. Hence ra*^^ ^ 5 x 10^'^ m~^ at Vq = 0. n* and as function of Vg can 
be modeled with, 

nl= / Die)fie,X\fil)de- D(6)/(e, 7^, /)cie 
Jo Jo 

nl = f D{e)[l-f{eX,f^l)]de- f D{e)[l - f{e,T^, fi')]de (8) 

J — oo J —oo 



and 



POO pO 

CBVG = e / Die)fie,%\fil)de-e / D ie)[l - fie,r,\ f^l)]de (9) 

Jo J-oo 

where / is the Fermi-Dirac distribution function, Te/h and /Ze/h are the respective carrier 
temperatures and Fermi levels. The superscript and 1 denotes the absence and presence 



of light excitation. D{e) = ^ + is the density-of-states, where eo is introduced 

to account for the electron-hold puddles. Due to the photo-excitation, the carriers will be 
driven away from equilibrium, characterized by a non-equilibrium Fermi energy /i^^^ and 
an elevated carrier temperatures 7^^^ compared to the ambient 7^°^. 
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At steady state, electrons and holes are allowed to thermalize among themselves, i.e. 
Te = Th and /ig = yujj, facilitated by femtosecond time scale carrier-carrier scattering 
processes [U E]. Here, T^j^ can be described by 7^^^ — % = QL/kq where Q is the absorbed 
laser power, L the device length and kq is the electronic thermal conductivity. Since a 
and kq are related through the Wiedemann-Franz relation, 7^^^ — To is then proportional 
to I/TqO", where the proportionality constant is determined to give us n*^^ ^ 5 x 10^^ m~^ 
at Vg = 0. This corresponds to T^)^ ~ % ~ 8K and 12 K at 7^ = 300 K and 200 K 
respectively. The photo-excited carriers n*^^ can then be numerically determined with Eq. 
j8jj9j by imposing charge conservation n* = n"^. Having calculated n*^i^ as a function of Vq 
then provides us with an estimatation of JpviVc) used in the main manuscript. 

In our analysis, we have extracted the photo-induced carrier density n* from electrical 
measurements described above. Alternatively, one can also estimates the photo-induced 
carrier density based on our light excitation condition. However, uncertainty in various 
parameters render it less accurate than the electrical method. Nevertheless, we can perform 
estimates of the photo-induced carrier density based on our light excitation condition. In 
our experiments, the laser power is P = 370 /iW with focal area aspot = f (0.7)^/im^. Light 
absorption at A=690nm (i.e. photon energy Eph = 1.8 eV) in graphene on 90nm Si02 is 
a ~ 2.5%. The photo-induced carrier density can be expressed as ra*^^ = MaPTrc/ Ephagpot, 
where M is the carrier multiplication factor and r^c is the carrier recombination time. Since 
^e/h ~ 5 X 10^'^ m~^, we estimate that Mr^c ~ 0.6 ps, which seems reasonable [B]. 



VI. INTRINSIC ELECTRON-PHONON LATTICE HEATING 

Electron-electron interaction results in an energy equilibration of the electronic system 
but does not lead to a net energy loss. The dominant energy loss pathways are due to 
phonons [3l [3-[9]. In particular, electronic cooling in graphene due to intrinsic acous- 
tic/optical phonon scattering processes has been well studied [TJ |8]. For example, the 
electron-lattice energy transfer mediated by acoustic phonons has the following power den- 
sity (Wm~^) given by [7], 




25 



where Dac ~ 20 eV is the acoustic phonon deformation potential and pm is mass density of 
graphene. For the experimental condition Te — Ti ~ 10 K and undoped graphene, Qac is only 
of the order of 10^ Wm~^. Under some doping and temperature conditions, the optical power 
density Qop may dominate over its acoustic counterpart however Q op /Qac is generally 
< 100 over the range of experimentally relevant conditions. 
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FIG. 1: (supplemental) Correction for the Photo Field-Effect, (a) Pliotocurrent amplitude and 

phase as a function of y-position (perpendicular to the graphene device) and gate voltage, (b) Fitting of the 

photo field-effect component of the photocurrent to Lorentzians for selected gate voltages. The gray-shaded 

area indicates the position of the f /im wide graphene device, which was excluded for fitting purposes, (c) 

Photocurrent as a function of y-position corrected for the photo field-effect for the same gate voltages as in 

(b). (d) Measured AC photocurrent (red) in the center of the graphene FET in Fig. 2 of the main text, and 

corresponding DC current (blue) as a function of gate voltage, (e) Photo field-effect at the center of the 

graphene channel (red) extracted from fits similar to the ones in (b). The photo field-effect is proportional 

to the transconductance (blue), (f) Photocurrent (red) corrected for the photo field-effect. The DC current 
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(blue) is plotted again as a reference. 
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FIG. 2: (supplemental) Spatial behavior of the AC photocurrent at low and high drain bias. 

Amplitude (a) and phase (b) of the AC photocurrent as a function of gate voltage and x-position along the 
graphene device for drain voltages from Vd^OV to -3V. 
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